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An  algorithm  for  finding  a  minimal  edge  coloring  of  a  bipartite  graph  in 
time  0(E  log  V)  is  presented*  Polynomial  time  algorithms  for  this  problem  have 
previously  been  given  by  Gabow  in  [1]  and  by  Gabow  and  Kariv  in  [2].  the  best 
time  bounds  being  0(E  log^  V)  and  0(V^  log  V). 

The  algorithm  is  based  on  using  fast  methods  for  finding  maximal 
matchings  in  semircgular  bipartite  graphs;  an  algorithm  for  finding  a  maximal 
matching  in  a  general  bipartite  graph  was  given  by  Hopcroft  and  Karp  in  [3]* 
Two  algorithms  for  finding  such  a  matching  are  given.  Although  the  second  one 
always  has  a  faster  running  time  of  0(max{E*  V  log  V  log^  D})»  the  first  one 
is  presented  for  the  sake  of  clarity. 

2.  KQimOH  ME  defihitions 

Throughout  this  paper  G  =  (V.E)  denotes  a  graph.  V  its  vertex  set  and  E 
its  edge  set.  G  =  (Vj.l^.E)  denotes  a  bipartite  graph  with  Vj  and  V2  being 
disjoint  vertex  sets  and  E  c  Vj  *  V2  being  the  edge  set.  D  denotes  the  maximal 
degree  of  any  vertex  in  V  =  Vj  u  V2. 

A  graph  is  said  to  be  regular  if  all  its  vertices  have  the  same  degree.  A 
bipartite  graph  is  said  to  be  semiregular  if  all  the  vertices  in  Vj  have  the 
same  degree  D.  the  maximal  degree  of  any  vertex  in  G;  it  is  said  to  be  higu- 
low  if  there  exists  an  integer  k  such  that  deg(v)  i  k  if  v  «  Vj  and  deg(v)  S  k 
if  v  e  Vj. 

An  Euler  partition  is  a  partition  of  the  edges  into  open  and  closed 
paths,  so  that  each  vertex  of  odd  degree  is  at  the  end  of  one  open  path*  and 
each  vertex  of  even  degree  is  at  the  end  of  no  open  paths. 

An  Euler  split  of  a  graph  G  =  (Vj,V2.E)  is  a  pair  of  graphs  Gj  = 
(Vj.l^.Ej)  and  G2  s  (V^.V2.E2)  where  Ej  and  E2  are  formed  from  an  Euler 
partition  of  E  by  placing  alternate  edges  of  each  path  into  Ej  and  E2 
respectively.  Any  vertex  of  even  degree  in  G  will  have  the  same  degree  in  both 
Gj  and  G2»  while  any  vertex  of  odd  degree  in  G  will  have  degrees  in  Gj  and  G2 
differing  by  one.  This  implies  that  if  G  is  semiregular*  and  D.  the  maximal 
degree  of  any  vertex  in  G  is  even,  then  both  Gj  and  G2  are  semiregular.  An 
algorithm  for  finding  an  Euler  split  in  time  0(E)  is  given  in  [1]. 

1.  ALGORITHM  1 


An  0(E  log  V)  time  algorithm  for  finding  a  maximal  matching  in  a 
semircgular  bipartite  graph  is  given.  The  algorithm  works  by  partitioning  E 
into  sets  Ej  and  E2  such  that  Gj  =  (Vj.V2.Ej)  and  G2  =  (Vj.V2«E2)  are  both 
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nontrivial  seniregular  bipartite  graphs.  If  the  graph  with  smaller  edge  set 
has  maximum  degree  one  it  is  the  required  matching;  otherwise  the  algorithm  is 
applied  recursively  to  that  graph.  Since  each  iteration  reduces  E  by  at  least 
a  factor  of  two.  the  algorithm  eventually  terminates. 

The  partitioning  procedure  to  obtain  the  semiregular  graphs  Cj  and  Gj  is 
as  follows.  An  Euler  split  of  G  is  made  giving  graphs  Gj  and  G2*  If'D,  the 
maximal  degree  of  any  vertex  in  G  is  even,  then  both  Gj  and  G2  are 
semiregular.  Otherwise  edges  are  moved  between  Gj  and  Gj  to  make  them 
seniregular.  This  is  described  more  precisely  below. 

Let  M  be  the  set  of  maximum  degree  vertices  in  G.  At  least  half  of  the 
vertices  in  M  will  have  even  degree  in  one  of  Gj  or  G2.  Without  loss  of 
generality  let  G2  be  that  graph  in  which  at  least  half  the  vertices  are  of 
even  degree.  Then  let  Mj  be  those  vertices  of  M  that  have  even  degree  in  G2 
and  M2  be  the  remaining  vertices  of  M. 

Next  an  Euler  split  of  G2  is  made  giving  graphs  Gjj  and  Gjj*  The 
vertices  of  Mj  have  the  same  degree  in  Gjj  and  (>22*  while  some  of  the  vertices 
in  Mj  have  even  degree  in  Gjj  and  odd  degree  in  G22  and  the  others  have  odd 
degree  in  Gjj  and  even  in  G22*  these  degrees  differ  by  one.  and  one  of  them  is 
the  degree  of  the  vertices  of  Mj  in  G2j  (and  in  ^27}*  Without  loss  of 
generality  let  G22  be  the  graph  in  which  at  least  half  the  vertices  of  M2  have 
even  degree.  Let  M21  be  the  subset  of  vertices  of  M2  that  have  even  degree  in 
Gj2  and  let  M22  be  the  remaining  vertices  of  H2« 

Now  one  of  the  graphs  Gjj  or  G22  is  combined  with  Gj  in  such  a  way  that 
vertices  in  Mj  will  have  the  same  degree  as  vertices  in  Kjj  in  the  combined 
graph.  The  new  graphs  are  named  Gj  and  G2  in  such  a  way  that  vertices  in  K22 
are  of  odd  degree  in  G2.  M^  and  M2  are  redefined  with  M2  reduced  in  size  by 
at  least  a  factor  of  two.  The  process  is  repeated  until  =  M  when  the 
vertices  in  M  all  have  the  same  even  degree  in  Gj.  The  partitioning  procedure 
is  shown  in  algol  like  form  below. 


PROCEDURE  PMHXIQfl 


G  =  (V.E) 

H  =  set  of  maximum  degree  vertices  of  G 
BEGIN 

Let  Gj»  G2  be  an  Euler  split  of  G; 

At  least  half  of  the  vertices  in  M  have  even  degree  in  one  of 
G1  or  G2.  Let  it  be  G2; 

Let  =  {v |v  €  H»  and  v  has  even  degree  in  G2); 

Let  M2  =  M  -  Mj8 
WHILE  (|2! 1 2  *  0)  DO 
BEGIN 

Let  G21,  G22  be  an  Euler  split  of  G22; 

Again  at  least  half  the  vertices  in  M2  have  even  degree  in 
either  G2j  or  G22*  Let  it  be  G22; 

Let  M2j  =  {v|v  €  M2»  and  v  has  even  degree  in  G22); 


Let  M22  ~  M2 


M 


21' 


IF  degree  of  a  vertex  in  Mj  in  G22 
then 


is  even 


else 


Gj  s=  Gj  u  C21.  G2  :=  G 


22 


END 


G2  :=  Gx  u  G22,  Gj  :=  G21; 
H|  #*  Mj  u  M2j*  M2  •*  M22; 


END 


CORRECTNESS 

It  is  necessary  to  show  that  all  the  vertices  in  M«  have  the  same  degree 
in  Gj  at  any  given  stage  of  the  algorithm*  and  likewise  in  G2*  The  same  result 
should  be  proven  for  vertices  in  1>2*  It  will  first  of  all  be  illustrated  by  an 
example* 

Consider  the  example  in  which  vertices  in  Mj  have  degree  5  in  G^  and 
degree  12  in  G2»  while  those  in  K2  have  degree  4  in  Gj  and  degree  13  in  G2* 

Then  vertices  in  have  degree  6  in  both  G2j  and  G22;  vertices  in 
have  degree  7  in  G^  and  degree  6  in  G22*  and  vertices  in  1I22  have  degree  6  in 
G2j  and  degree  7  in  G22*  So  the  assignments  Gj  =  Gj  u  G2j»  G2  =  G22*  M*  =  M| 
u  M2j*  and  M2  *  M22  are  made* 

Now  vertices  in  Mj  have  degree  11  in  Gj  and  degree  6  in  G2*  and  vertices 
in  M2  have  degree  10  in  Gj  and  degree  7  in  G2* 

By  considering  respectively  the  cases  in  which  the  degree  in  Gj  of 


vertices  in  Mj  is  one  greater  or  one  lesser  than  that  of  vertices  from  Hj  it 
can  be  proven  by  induction  that  the  degree  of  vertices  in  Mj  is  the  sane  in 
each  of  Gj  and  G2  and  likewise  for  Mj.  Thus  all  the  vertices  in  M  have  the 
same  degree  in  each  of  Gj  and  G2  when  the  partitioning  procedure  terminates. 

To  show  that  Gj  and  G2  are  both  semiregular  it  is  necessary  to  show  that: 

deg(v)  in  Gj,  v  not  in  K  £  deg(v)  in  Gj»  v  e  M»  and  similarly  in  Gj. 

This  is  proven  by  an  induction  using  the  inductive  hypothesis  that: 

deg(v)  in  Gj,  v  not  in  M  5  deg(v)  in  Gj.  v  €  Hj,  and  likewise  in  Gj. 

To  show  that  both  G^  and  G2  are  nontrivial  the  following  inductive 
hypothesis  is  proven: 

deg  of  vertices  of  Mj  in  G^  >  0. 

In  fact  this  degree  is  even  and  so  the  degree  of  vertices  from  Mj  in  G^j  and 
G22  is  nonzero*  The  induction  now  follows. 

TIHIHG 

Each  iteration  of  the  while  loop  reduces  the  size  of  Mj  by  at  least  a 
factor  of  two.  So  after  at  most  0(log  V^)  iterations  iMjl  -  0  and  the 
procedure  terminates.  Each  iteration  of  the  while  loop  takes  time  0(E).  So  to 
obtain  G^  and  G2  takes  time  0(E  log  Vj)  £  0(E  log  E/D).  Thus  the  time  T(E) 
taken  to  find  a  matching  is  given  by: 

T(E)  =  0(E  log  E/D)  +  T(E/2)  =  0(E  log  E/D)  :S  0(E  log  V). 

In  fact  this  algorithm  will  produce  a  matching  covering  all  the  vertices 
of  maximal  degree  in  a  general  bipartite  graph  in  time  0(E  log  E/D). 

4.  ALGORITHM  2 

An  0(E  ♦  V  log  V  log  D)  time  algorithm  for  finding  a  maximal  matching  in 
a  semiregular  bipartite  graph  is  given. 

OBSEmnON 

It  is  known  that  if  a  network  has  a  maximal  flow*  with  the  flow  through 
each  vertex  being  integral,  then  it  has  a  maximal  flow  such  that  the  flow 
through  each  edge  is  integral  [4.  pll3). 


In  particular*  if  the  edges  of  a  bipartite  graph  are  assigned  positive 
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weights*  so  that  the  sum  of  the  weights  at  each  vertex  is  at  most  one*  and  the 
weight  at  each  vertex  of  is  one*  then  the  graph  has  a  matching  covering 
every  vertex  of  Vj. 

By  shifting  the  weights  between  edges*  while  maintaining  a  constant 
weight  at  each  vertex*  and  deleting  edges  of  weight  zero*  a  new  smaller  graph 
is  obtained  containing  just  as  large  a  matching. 

METHOD 

Vertices  of  small  degree  are  merged  together  so  that  all  vertices  have 
degrees  between  D/2  and  D.  Dow  V*D  s  0(E).  This  simplifies  the  timing 
analysis.  Every  edge  is  given  weight  1/D.  Using  depth  first  search*  cycles  in 
the  graph  among  edges  of  weight  1/D  are  found.  When  a  cycle  is  found  alternate 
edges  in  the  cycle  are  deleted;  the  other  edges  in  the  cycle  have  their  weight 
doubled.  This  is  continued  until  there  are  no  cycles  among  edges  of  weight 
1/D.  Cycles  are  then  found  among  edges  of  weight  2/D*  A/D  ..in  turn  until  no 
further  increase  in  edge  weights  can  be  obtained  in  this  way. 

A  graph  with  at  most  0(V  log  E/V)  weighted  edges  is  obtained*  such  that 
the  sum  of  the  weights  at  each  vertex  in  is  one.  Algorithm  1  is  now 
adapted  to  find  this  matching. 

Each  edge  is  considered  to  have  multiplicity  D  times  its  weight  .  Four 
copies  of  each  edge  are  kept*  one  in  each  of  G^*  62*  ^21  an<*  ^22*  ^ken  asking 
an  Euler  split*  each  edge  added  to  an  Euler  path  is  made  to  occur  as  often  as 
it  can  at  that  point  in  the  path;  the  multiplicities  of  the  copies  of  the  edge 
are  changed  accordingly.  Otherwise  one  proceeds  just  as  in  algorithm  1. 

As  with  algorithm  1  this  algorithm  can  be  used  to  find  a  matching 
covering  all  the  vertices  of  maximal  degree  in  a  general  bipartite  graph  in 
the  same  time  bound. 

TTMIHC 

One  iteration  of  the  procedure  in  algorithm  1  cuts  the  number  of  edges  in 
half  (counting  edges  according  to  their  multiplicity)*  but  may  well  not  reduce 
the  size  of  the  structure  stored*  affecting  the  timing  given  with  algorithm  1. 

To  find  the  graph  of  size  0(V  log  E/V)*  0(E)  time  is  needed.  To  use 
algorithm  1  one  requires  time  0(V  log  E/V  log  E/D)  to  halve  the  number  of 
edges  (counted  according  to  their  multiplicity)*  and  time  0(V  log  E/V  log  E/D 

O 

log  D)  =  0(V  log  V  log  E/V)  to  obtain  a  maximal  matching. 

2 

Thus  the  overall  time  taken  is  0(E)  for  E  £  0(Vlog  V  (loglog  V))  and 


0(V  log  V  log2  E/V)  otherwise*  which  is  always  better  than  the  0(E  log  V)  of 
algorithm  1. 

1.  COLORING  HIE  EDGES  A  BIPARTITE  £M£U. 

An  0(E  log  V)  tine  algorithm  for  finding  a  minimal  edge  coloring  is 
given*  It  is  minimal  in  the  sense  that  the  fewest  possible  colors  are  used.  It 
is  shown  in  Cl]  that  this  is  D  colors. 

The  algorithm  is  based  on  divide  and  conquer.  When  a  graph  has  vertices 
of  even  maximal  degree*  using  the  method  of  Euler  partitions  it  is  split  into 
two  subgraphs  of  equal  maximal  degree  which  are  then  recursively  colored.  On 
occasion  when  the  graph  has  vertices  of  odd  maximal  degree  a  matching  M 

covering  all  the  vertices  of  maximal  degree  has  to  be  found.  This  is  obtained 

by  using  algorithm  2. 

ALGQ&IIflll 

If  D  is  odd  find  a  matching  as  described  above;  color  it  and  delete  it 
from  G.  Set  D  =  D  -  1 . 

Hake  an  Euler  split  of  G  to  give  two  bipartite  graphs  Gj  and  G2  each 

having  vertices  of  maximal  degree  D/2. 

WLOG  assume  G,  has  a  smaller  edge  set  than  G,  (otherwise  swap  the  labels 
1  k  ^  k+1 

Gj  and  62)*  and  recursively  color  Gj*  Let  2  <  D/2  s  2  -  r*  Add  the  r 

smallest  sets  of  colored  edges  to  G2»  delete  them  from  the  set  of  colored 

edges*  and  recursively  color  G2* 

A  similar  method  was  used  in  [3]  and  led  to  the  current  presentation  of 
the  algorithm*  That  exactly  D  colors  are  used  can  be  shown  by  induction* 

TIMING 

Excluding  the  time  taken  to  find  the  matchings*  the  time  taken  is  given 

by: 

T(E.  D)  =  T(Ej,  ID/2J)  +  T(E2  u  E3*  LD/2J  +  r)  +  0(E)* 

where  Ej  u  Ej  s  E  and  E3  is  the  union  of  the  r  sets  of  colored  edges  added  to 
G2.  For  D  a  power  of  two  T(E,  D)  =  o(E  log  D).  In  all  other  cases  as  |£3 1  £ 
2r  lEjI/D  £  I E2 I  one  finds  that  T(E*  D)  =  0(2E  log(D  +  r))  =  0(E  log  D). 

The  time  required  for  finding  the  matchings  is  bounded  by  0(max(E  log  D* 
V  logV  log  D})  £  0(max{E  log  D*  E  log  V})  and  hence  the  total  time  required 
is  bounded  by  0(E  maxllog  D.  log  V))  which  is  0(E  log  V).  For  a  graph  without 


multiple  edges  a  time  bound  of  0(E  max {log  D*  log  V})  is  obtained 


L.  MATCHINGS  Ill  HIGH  L Oil  GRAPHS 

By  pruning  the  high-low  graph  a  semiregular  graph  with  D  =  k  can  be 
obtained.  D  being  the  maximal  degree  of  any  vertex  in  the  semiregular  graph* 
The  matching  algorithm  is  then  applied  to  this  graph  to  obtain  a  maximal 
matching  for  the  high-low  graph*  So  a  maximal  matching  in  a  high-low  graph  can 
be  found  in  time  0(max{E*  V  log  V  log2  E/V})*  High-low  graphs  were  defined  in 
[3]  and  the  above  method  for  pruning  was  described  there* 
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